If we are only using Paracou Biodiversity and Control plots we lack of Symphonia globulifera morphotype to have a balanced stratified sampling (\(\frac{99}{400}\) instead of \(\frac{99}{400}\), see Table 1.1), and we only have 6 S. indet which means we can’t hope they are in fact S. globulifera. Including treatment plots only allow us to reach 168 S. globulifera morphotype but with a lot of S. indet in which we might find other S. globulifera and thus reached a balanced sampling. We observed a strong relation to be tested between water table depth and S. globulifera presence (see Figure 1.1). So we may estimate with the layer the number of S. globulifera alive in Paracou. But including treatment plots necessitate to control for the light variable due to anthropic gaps of the logging. We thus need to test the effect of the treatments and gaps on individuals growth, mortality, and recruitment to evaluate the effect of the treatment on Symphonia Paracou population. So in parrallel we are looking at other Guyafor plots looking for Symphonia globulifera with both diameters inventories for at least three censuses and environmental data.
| espece | B | C | T1 | T2 | T3 | Total |
|---|---|---|---|---|---|---|
| globulifera | 65 | 34 | 20 | 31 | 18 | 168 |
| Indet. | 6 | 165 | 166 | 66 | 0 | 403 |
| sp.1 | 367 | 62 | 81 | 102 | 127 | 739 |
Table 2.1 shows few individuals in other Guyafor plots. It is due to the fact that most of Guyafor network plots include few bottom lands (Pascal Petronelli, personnal communication). Consequently it seems if we want both good quality diameter inventories and environment data we will have to sample all individuals in Paracou. The question stay to sample or not in treatment plots ? And more generally the question, is about to know if we will have sufficient environmental data (gaps map, lidar) to control anthropic and natural gaps effects ?
| NomForet | globulifera | Indet. | sp.1 |
|---|---|---|---|
| Paracou | 252 | 666 | 868 |
| Saül Diadema (Limonade) | 17 | 0 | 0 |
| Nouragues | 14 | 1 | 16 |
| Acarouany (Javouhey) | 1 | 1 | 12 |
| Régina St Georges | 1 | 31 | 6 |
We gathered genetic material (ddRADseq) of french Guiana from Torroba-Balmori unpublished data (Paracou and Regina). We cleaned
fastqfiles after a check withfastQCheckallowing us to correct two sequences by removing theim for individuals \(PR_{49}\) and \(RG_1\). We usedipyradfor the interactive assembly of ddRADseq data sets on genotoul cluster (with denovo assembly, AATT and AT restriction overhang, 85% clustering threshold and a minimum of 48 sample per locus).
We used
vcfRto load SNPs data into R, and we transform it in genligh object foradegenet. We related indivdual IDs to their population and coordinates with links table. We coded population in 4 subset for Symphonia globulifera and sp1 in both Paracou and Régina (\(PR_{gl}\), \(PR_{sp}\), \(RG_{gl}\), \(RG_{sp}\)). Population definition was used to transform vcf file to structure file with PGDspider for further genetic structure analysis with STRUCTURE software. We corrected and transformed in UTM coordinates to compute kinship distance matrix with SPAGEDI.
For the moment habitat association are tested only with Water Table Depth but should further be tested for different environmental variable. We can see that S. globulifera is significantly positively associated to water table depth under 60 cm on most of cases (11 plots over 15 for morphotype and 2 plots for genotypes, knowing that the 4 lacking plots for morphotype are certainly due to too coarse rasterization of the habitat shapefile). And S. globulifera is also significantly negatively associated to water table depth over 100 cm (12 plots over 15 for morphotypes and 2 plots for genotypes). Last but not least S. sp.1 does not show significant association for genotypes and only 3 for morphotype (2 positive association to water table depth over 100 cm and one negative association to water table depth under 60 cm). Finally, absence of a lot of significant association for genotypes might not be due to an absence of association but more to a lack of power. Indeed only individuals associated more than 90% to a genotype in STRUCTURE analysis were kept resulting in relatively few individuals.
| WTD | species | plot | pvalue | two.tail.test |
|---|---|---|---|---|
| Nappe au dela de 100 cm | globulifera | 1 | 0.0027624 | - ** |
| Nappe au dela de 100 cm | globulifera | 2 | 0.0027624 | - ** |
| Nappe au dela de 100 cm | globulifera | 4 | 0.0027624 | - ** |
| Nappe au dela de 100 cm | globulifera | 6 | 0.0027624 | - ** |
| Nappe au dela de 100 cm | globulifera | 7 | 0.0027624 | - ** |
| Nappe au dela de 100 cm | globulifera | 8 | 0.0027624 | - ** |
| Nappe au dela de 100 cm | globulifera | 11 | 0.0027624 | - ** |
| Nappe au dela de 100 cm | globulifera | 13 | 0.0027624 | - ** |
| Nappe au dela de 100 cm | globulifera | 14 | 0.0027624 | - ** |
| Nappe au dela de 100 cm | globulifera | 15 | 0.0027624 | - ** |
| Nappe entre 0 et 60 cm | globulifera | 1 | 1.0000000 | + *** |
| Nappe entre 0 et 60 cm | globulifera | 3 | 1.0000000 | + *** |
| Nappe entre 0 et 60 cm | globulifera | 4 | 1.0000000 | + *** |
| Nappe entre 0 et 60 cm | globulifera | 5 | 1.0000000 | + *** |
| Nappe entre 0 et 60 cm | globulifera | 7 | 1.0000000 | + *** |
| Nappe entre 0 et 60 cm | globulifera | 8 | 1.0000000 | + *** |
| Nappe entre 0 et 60 cm | globulifera | 10 | 1.0000000 | + *** |
| Nappe entre 0 et 60 cm | globulifera | 11 | 1.0000000 | + *** |
| Nappe entre 0 et 60 cm | globulifera | 13 | 1.0000000 | + *** |
| Nappe entre 0 et 60 cm | globulifera | 14 | 1.0000000 | + *** |
| Nappe entre 0 et 60 cm | globulifera | 15 | 1.0000000 | + *** |
| Nappe entre 60 et 100 cm | globulifera | 6 | 1.0000000 | + *** |
| Nappe entre 60 et 100 cm | globulifera | 13 | 1.0000000 | + *** |
| Nappe au dela de 100 cm | sp.1 | 3 | 1.0000000 | + *** |
| Nappe au dela de 100 cm | sp.1 | 15 | 1.0000000 | + *** |
| Nappe entre 0 et 60 cm | sp.1 | 1 | 0.0027624 | - ** |
| WTD | species | plot | pvalue | two.tail.test |
|---|---|---|---|---|
| Nappe au dela de 100 cm | globulifera | 4 | 0.0027624 | - ** |
| Nappe au dela de 100 cm | globulifera | 5 | 0.0027624 | - ** |
| Nappe entre 0 et 60 cm | globulifera | 4 | 1.0000000 | + *** |
| Nappe entre 0 et 60 cm | globulifera | 5 | 1.0000000 | + *** |
We still do not have BRIDGE raw data too compare intraindividual versus interindividual functional trait variations.
I wanted to test for an eventual effect of logging through light and disturbance on Symphonia individuals growth. I used control and treatment plots between 1988 and 1992. I kept only trees already present in 1988 and still alive in 1992, and calculate their growth during this time. I then tried to look at the effect of both distance to the closest logging gaps (\(d_{gaps}\) in \(m\)), original dbh of inidividuals in 1998 (\(dbh_{1998}\) in \(cm\)), and their interaction on growth (\(growth\) in \(cm\)). I used a bayesian model following a log normal law for growth : \[growth \sim log \mathcal{N}(\alpha*log(d_{gaps}+1) + \beta*dbh_{1998} + \gamma*log(d_{gaps}+1)*dbh_{1998}, \sigma) \] The model seems to have correctly converged besides log likelyhood seems to be constrained to a maximal value (see markov chain plot), and parameters seems not much correlated besides a small link between \(\beta\) and \(\gamma\) (see parameters pairs plot). So if we consider the model as valid, parameter posterior distribution seems to indicate a strong effect for the logarithm distance to the gap (\(\alpha\) distribution don’t overlap 0 with a mean at \(\alpha_m = 0.18\)). The growth of Symphonia individuals being increased close to gaps by almost 14% (\(e^{\alpha_m}=1.14\)).
\[growth = e^{0.13*log(d_{gaps}+1) + 0.03*dbh_{1988} + -0.01*log(d_{gaps}+1)*dbh_{1988}} \]